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Abstract — Typically,  actuator  dynamics  are  ignored  when  de¬ 
signing  flight  control  allocators  for  aircraft  because  the  band- 
widths  of  actuators  are  normally  much  larger  than  the  fre¬ 
quencies  of  the  vehicle’s  rigid  body  modes.  Unfortunately, 
this  in  not  always  the  case,  particularly  when  dealing  with 
non-aerodynamic  surface  actuators.  Ignoring  the  interactions 
between  constrained  control  allocators  and  actuator  dynamics 
can  have  serious  consequences.  In  this  work,  a  method,  which 
post-processes  the  output  of  a  control  allocation  algorithm,  is 
developed  to  compensate  for  actuator  dynamics.  The  actua¬ 
tors  can  have  dynamics  which  are  either  first-order,  second- 
order  with  no  zeros,  or  second-order  with  a  single  zero.  The 
method  developed  here  solves  for  a  gain,  which  multiplies 
the  commanded  change  in  control  effector  setting  as  com¬ 
puted  by  the  control  allocator.  This  approach  is  not  computa¬ 
tionally  intensive  and  thus  has  the  added  benefit  of  being  an 
algorithm  which  can  operate  in  real-time  on  a  typical  flight 
computer.  Likewise,  this  approach  is  applicable  to  both  the 
saturated  and  unsaturated  control  effector  cases.  The  basic 
premise  of  this  method  is  to  post-process  the  output  of  the 
control  allocation  algorithm  to  overdrive  the  actuators  so  that 
at  the  end  of  a  sampling  interval,  the  actual  actuator  positions 
are  equivalent  to  the  desired  actuator  positions. 
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1.  Introduction 

Control  allocation  algorithms  are  used  to  compute  aircraft 
control  surface  deflections  which  produce  a  desired  set  of  mo¬ 
ment  or  acceleration  commands.  A  number  of  approaches 
have  been  developed  that  ensure  that  the  commands  provided 
to  the  effectors  are  physically  realizable.  These  actuator  com¬ 
mand  signals  are  feasible  in  the  sense  that  they  do  not  exceed 
hardware  rate  and  position  limits.  Buffington[l]  and  Bod- 
son  [2]  developed  methods  which  take  into  account  rate  and 
position  limits. 

A  review  of  the  constrained  control  allocation  literature 
shows  that  the  coupling  effects  that  result  from  combin¬ 
ing  constrained  control  allocators  and  actuator  dynamics  has 
largely  been  ignored,  although  some  recent  research  has  been 
performed  on  directly  including  actuator  dynamics  in  the  con¬ 
trol  allocation  problem[3],  [4].  The  underlying  assumption  of 
most  previous  work  is  that  actuators  respond  instantaneously 
to  commands.  This  assumption  may  at  first  seem  justified  be¬ 
cause  in  practice,  actuator  dynamics  are  typically  much  faster 
than  the  rigid  body  modes  that  are  to  be  controlled.  How¬ 
ever,  interactions  between  a  constrained  control  allocator  and 
an  actuator  with  linear  dynamics  can  result  in  a  system  that 
falls  well  short  of  its  potential.  As  an  example,  a  simula¬ 
tion  was  run  with  a  linear  programming  based  control  allo¬ 
cation  algorithm  mixing  four  control  effectors  to  obtain  a  de¬ 
sired  set  of  moments,  d des  £  M3.  In  this  simulation  run, 
the  actuator  dynamics  for  each  control  effector  were  set  to 
§  =  ^5  and  the  commanded  effector  positions,  as 

computed  by  the  control  allocator,  were  used  to  initialize  the 
allocator  at  the  next  timestep.  To  illustrate  the  effects  of  con¬ 
strained  control  allocator  and  actuator  dynamics  interactions, 
consider  Figures  1  and  2.  The  objective  is  to  make  the  accel¬ 
erations  produced  by  the  control  effectors  ( B8 )  equal  to  the 
commanded  accelerations  (d^es).  Figure  1  shows  that  when 
there  are  no  actuator  dynamics,  the  desired  result  is  achieved, 
namely  B8  =  d^es.  When  actuator  dynamics  are  included, 
the  results  are  as  shown  in  Figure  2  where  it  is  obvious  that 
B8  /  d des-  The  goal  of  this  work  is  to  develop  a  scheme, 
which  compensates  for  the  effects  of  actuator  dynamics,  so 
that  even  when  actuator  dynamics  are  included,  results  like 
those  shown  in  Figure  1  are  still  obtained.  As  was  shown 
in  the  example  above,  the  control  allocator/actuator  interac- 
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Figure  1.  Acceleration  Commands  and  Accelerations 
Produced  by  the  Controls  -  No  Actuator  Dynamics. 
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Figure  3.  Control  allocator  and  actuator  interconnection. 
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Figure  2.  Acceleration  Commands  and  Accelerations 
Produced  by  the  Controls  -  Actuator  Dynamics  Present. 


tion  can  yield  control  effector  positions  which  produce  signif¬ 
icantly  different  accelerations  than  the  commanded  accelera¬ 
tions.  Bolling [5]  has  shown  that  the  interaction  between  first- 
order  actuator  dynamics  and  constrained  control  allocation 
algorithms  can  be  eliminated  by  overdriving  the  actuators.  In 
previous  work,  [6]  the  effects  of  the  control  allocator/actuator 
interaction  for  first-order  and  second-order  (with  no  zeros) 
actuator  dynamics  has  been  taken  into  account.  That  work 
is  expanded  here  to  include  second-order  actuator  dynamics 
with  a  single  zero.  Having  the  first-order,  second-order,  and 
second-order  with  a  zero  cases  defined,  this  work  makes  it 
possible  to  create  higher-order  actuators  as  a  combination  of 
these  simpler  forms. 


tor  dynamics  on  control  system.  Following  this,  the  details 
of  the  interaction  between  constrained  control  allocators  and 
second-order  actuator  dynamics  with  a  zero  are  provided  in 
detail.  Simulation  results  are  presented  for  a  four  control  sur¬ 
face  vehicle  using  a  linear  programming  based  control  allo¬ 
cation  algorithm  that  takes  into  account  control  effector  posi¬ 
tion  and  rate  limits.  The  results  show  that  the  desired  control 
effector  positions  can  be  obtained  by  post-processing  the  con¬ 
trol  allocation  commands. 

2.  System  Definition 

Figure  3  shows  the  system  to  be  analyzed  in  this  work.  Inputs 
to  the  control  allocation  algorithm  consist  of  a  vector  of  de¬ 
sired  moment  or  acceleration  commands,  d des  £  Mn,  and 
a  vector  containing  the  current  control  surface  deflections, 
S  £  Mm.  The  output  of  the  control  allocator  is  the  com¬ 
manded  control  surface  deflection  vector,  Scrnd  £  Mm.  The 
actuator  dynamics  respond  to  Scmci  to  produce  the  actual  con¬ 
trol  deflections,  S.  The  individual  actuators  are  assumed  to 
have  uncoupled  dynamics,  hardware  rate  limits  ±£max,  and 
position  limits  Sm in,  <5max.  In  most  control  allocator  imple¬ 
mentations,  rate  limits  are  taken  into  account  by  converting 
them  into  effective  position  limits  at  the  end  of  the  next  sam¬ 
pling  period  and  constraining  the  effector  commands  to  re¬ 
spect  the  most  restrictive  of  the  rate  or  position  limits,  i.e., 

S  =  min(<5max,  S  -f  ^max^)  /^\ 

(5  max(<5min,  S  ^maxAf) 

where  d  is  the  current  location  of  the  control  effectors,  d 
are  the  most  restrictive  upper  and  lower  bounds  on  the  effec¬ 
tors,  respectively,  and  At  is  the  sampling  period  of  the  digital 
flight  control  system. 

A  few  comments  are  in  order  regarding  Figure  3  and  the  anal¬ 
ysis  described  above.  First,  note  that  the  instantaneous  posi¬ 
tion  limit  given  by  Equation  1  makes  use  of  a  sampled  vector 
of  actuator  position  measurements  to  compute  the  maximum 


In  this  work,  results  for  the  first-order  and  simple  second- 
order  cases  are  provided  to  illustrate  the  effects  of  actua¬ 
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distance  that  the  actuator  can  move  in  the  next  time  instant. 
This  is  in  contrast  to  using  the  previous  value  of  actuator  com¬ 
mand  vector,  5Cmd ,  as  is  often  done  in  simulation.  The  mo¬ 
tivation  for  using  actuator  measurements  is  that  when  actua¬ 
tor  dynamics,  disturbances,  and  uncertainties  are  taken  into 
account,  the  actuator  command  vector  and  the  true  actuator 
positions  will  differ.  This  difference  can  cause  a  control  allo¬ 
cator  to  generate  inappropriate  actuator  commands  that  do  not 
deliver  the  desired  moments  or  accelerations.  Thus,  feeding 
the  measured  actuator  position  vector  to  the  control  allocator 
has  the  advantage  of  reducing  uncertainty  in  the  actuator  po¬ 
sition.  It  is  important;  however,  to  be  aware  of  the  effects  of 
using  the  measured  actuator  position  vector. 

3.  Attenuation  of  Zero-Order-Hold 
Inputs  For  First-Order  and 
Second-Order  (No  Zeros)  Actuator 
Dynamics 

Referring  to  Figure  3,  the  desired  situation  would  be  for 
S  =  Scrnd .  However,  actuator  dynamics  alter  the  command 
signals  so  that,  in  general,  S  ^  Scmd •  For  actuators  with  high 
bandwidths  relative  to  the  rigid  body  modes,  this  is  not  a  se¬ 
rious  concern.  However,  situations  exist  where  the  actuator 
dynamics  are  not  sufficiently  fast  and  need  to  be  taken  into 
account.  In  this  section,  the  effects  of  first-order  actuator  dy¬ 
namics  and  simple  second-order  (no  zero)  actuator  dynamics 
on  the  system  shown  in  Figure  3  will  be  discussed.  Let  the 
dynamics  of  a  single  actuator  be  represented  by  a  continuous 
time  first-order  transfer  function  of  the  form 

M  =  (2) 

^cmd(^)  S  “h  CL 

The  discrete  time  representation  of  the  first-order  actuator  dy¬ 
namics  equation  is  given  by 

S(tk+i)  =  ®5(tk)  +  r5cmd(tk)  (3) 

where  4>  =  e~aAt,  T  =  1  —  e~aAt ,  and  it  has  been  assumed 
that  the  input  to  the  actuator  dynamics,  Scrnd(tk ),  is  held  con¬ 
stant  over  each  sampling  period.  The  command  to  the  actua¬ 
tor  can  be  written  as 

$cmd(pk)  =  A5crndCA  (Pk)  T~  ^(pk)  (4) 

where  the  commanded  incremental  change  in  actuator  po¬ 
sition  over  one  timestep  is  defined  by  A ScrndCA(tk)  = 
$ cmdc a  (tk)  -S(tk)  and  where  ScrndcA  (tk)  is  the  actuator  po¬ 
sition  command  from  the  control  allocator.  Since  the  effec¬ 
tor  commands  are  held  constant  over  each  sampling  period, 
A S cmdc A^k)  wiH  appear  to  the  actuators  to  be  a  step  com¬ 
mand  from  the  measured  position.  Recall  that  ScrndCA  (tk)  is 
calculated  by  the  control  allocator  based  upon  the  assumption 
that  the  effector  will  respond  instantaneously  to  commands. 
Substituting  Equation  4  into  Equation  3  yields 

3{tk+ 1)  =  mtk)  +  r  [A ScrndCA(tk)  +  S(tk)]  (5) 

Since  T  <  1,  the  incremental  command  signal  from  the  con¬ 
trol  allocation  algorithm,  A ScrndCA(tk)  is  attenuated  by  the 


Figure  4.  Block  diagram  of  command  increment 
compensation. 


actuator  dynamics,  thus  5(tk+i)  ^  ScmdCA(Pk)-  The  objec¬ 
tive  is  to  find  a  gain,  M,  that  modifies  the  output  of  the  con¬ 
trol  allocation  algorithm  such  that  5(tk+ i)  =  5 Crude A^Pk)  = 
crude A^Pk)  "F  fi(pk)*  Hence, 

5{tk+ i)  =  ®5(tk)  +  r  [MAS cmdc A  (tk)  +  S(tk)]  (6) 

and,  solving  for  M  yields 


(7) 


Thus  the  actuator  command  signal  must  be  modified  such  that 


3 cmdipk )  —  MA ScmdcA^Pk)  "F 

1  (8) 

—  p  A  $cmdc  A^Pk)  3“  $(pk) 

Replacing  Scmd(tk)  in  Equation  3  with  5crnd{tk)  from  Equa¬ 
tion  8  yields 


3(pk+l)  —  *&3(pk')  "F  r foemdipk) 


=  QS(tk)  +  E 


i 

r 


A ficmdcA^Pk)  T  3(pk) 


(9) 


which  yields  the  desired  actuator  position. 


Since  T  can  be  computed  from  the  known  quantities  a  and 
At,  one  can  compensate  for  command  increment  attenuation 
using  Equation  8.  For  a  bank  of  decoupled  first-order  actua¬ 
tors  with  nominal  bandwidths  of  a$,  corresponding  values  of 
Ti  can  be  computed  using  T^  =  (1  —  e~aiAt).  The  com¬ 
mand  increment  compensation  can  then  be  implemented  in 
discrete  time  as  shown  in  Figure  4.  In  Figure  4,  the  dashed 
lines  are  for  compensation  for  second-order  actuator  dynam¬ 
ics,  as  will  be  discussed  shortly,  and  therefore,  are  not  part  of 
the  first-order  actuator  dynamics  compensation  scheme.  Note 
that,  for  multiple  actuators,  M  in  Figure  4  is  a  diagonal  matrix 
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with  the  entries  on  the  main  diagonal  being  Ti,  T2,  •  •  •  ,  Tm, 
where  the  subscript  m  is  defined  as  the  number  of  control  ef¬ 
fectors.  Hence,  the  magnitude  of  the  control  allocation  com¬ 
mand  increment  is  modified  to  counteract  the  attenuation  that 
results  from  the  interaction  between  first-order  actuator  dy¬ 
namics  and  the  control  allocator. 

With  the  knowledge  gained  from  the  first-order  actuator  dy¬ 
namics  case,  it  is  now  possible  to  extend  the  previous  results 
to  second-order  actuator  dynamics.  Let  the  second-order  ac¬ 
tuator  dynamics  be  represented  by  the  general  equation 


(5(s) 
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Scmd(s)  S 2  +  2 (,U)nS  ■ 

which,  in  state-space  form,  becomes 
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Using  Equations  11  and  12,  the  discrete-time  solution  to 
the  second-order  actuator  dynamics  differential  equation  be¬ 
comes 


S(tk+ll  =  [  $1,1  $1,2 


s(tk) 

Hh) 


rtk+i 


■SCmd(tk)  2(4+i  -  r)dr 

Jtk 


where 


<j>  -  gA(tfe_)_i 


$1,1  $1,2 

$2,1  $2,2 


(13) 


(14) 


is  the  state  transition  matrix.  Since  A  is  constant,  the  state 
transition  matrix  depends  on  the  time  difference,  At  = 
tk+i  —  tk,  and  the  explicit  dependence  on  time  has  been  elim¬ 
inated.  Performing  the  operations  required  in  Equations  13 
and  14  produces 


S(tk+%)  —  Ci5(tk)  +  C2S(tk)  +  C^cmd(tk)  (15) 


where 


n  ~  n  , 

Ci  =  — e 


^  _crA t 


sm 


Ud 


uOdAt  +  arctan  ( 

—  (7 


7  At 
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Ud 
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odd  =  \/l  —  C2  an(l  a  =  —  C^n-  Substituting  Equation  4 

into  Equation  15  produces 

S(tk&)  =  C\8(tk)  +  C2S(tk)  +  Cs  [A ScmdCA(tk)  +  S(tk)] 

(17) 

The  objective  is  to  find  a  gain  M  that  will  modify 
A ScmdcA^k)  in  such  a  way  that  6(tk+ 1)  —  $ crude  a(^  k)- 
Hence,  it  is  desired  to  find  M  such  that 


A  ScmdcA^k)  +  S(tk)  —  Ci5(tk)  +  C25(tk)  + 

C3  [MA ScmdcA(tk)  +  5(tk)] 

Solving  for  M  gives 

A ScmdcA^k)  +  (1  ~  C3  —  Ci)S(tk)  ~  C2S(tk) 


(18) 


Mm 


C3A ScrndcA  (^fc) 


(19) 

Clearly,  unlike  the  simple  first-order  case,  this  compensation 
is  more  complex.  In  fact,  this  requires  not  only  the  position 
of  the  control  effector,  but  also  the  rate  of  change  of  the  con¬ 
trol  effector.  In  practice,  an  estimator  would  be  designed  (for 
example,  a  Kalman  filter)  to  estimate  5(tj. c).  Referring  to  Fig¬ 
ure  4,  for  a  system  with  more  than  one  actuator,  M  would  be 
a  diagonal  matrix  with  entries  along  the  main  diagonal  be¬ 
ing  Mi,  M2,  •  •  •  ,  Mm.  Here,  Mi  would  be  computed  using 
Equation  19  and  cjn. ,  Q  corresponding  to  the  ith  actuator. 

4.  Attenuation  of  Zero-Order-Hold 
Inputs  for  Second-Order  Actuator 
Dynamics  With  A  Single  Zero 

The  results  of  the  previous  section  are  now  extended  to  the 
case  of  second-order  actuator  dynamics  with  a  single  zero. 
Let  the  actuator  dynamics  be  represented  by 


S(s) 


k(s  +  a) 


3cmd(s)  32  +  2(^idnS  H~  L d2^ 

In  this  case,  it  is  easily  seen  that  Equation  20  contains  deriva¬ 
tives  of  the  input  signal.  In  order  to  eliminate  Scmd(t ),  intro¬ 
duce  an  intermediate  variable,  z(s),  so  that 


k(s  +  a) 


s2  +  2(cjns  +  u;2 


and  let 
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Then,  it  is  easily  seen  that  Equation  22,  in  state- space  form, 
becomes 
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(25) 

From  Equation  23,  the  output  equation  becomes 


S(t)  =  kaz(t)  +  kz(t)  =  ka  k 


Equations  24  and  25  describe  the  second-order  actuator  dy¬ 
namics. 

Using  Equations  24  and  25,  the  discrete-time  solution  to 
the  second-order  actuator  dynamics  differential  equation  be¬ 
comes 

6(4+ 1)  =  (ka$  1,1  +k$2,i)z(tk) 

+  (fca$ +2  +  fed? 2,2)  z(tk)-\- 

rtk+i 

5cmd(tk )  ka$i>2{tk+i  -  T)dr+  (26) 

Jtk 

ptk  +  1 

5cmd(tk)  2(4+1.- r)d,T 

Jtk 

where  was  defined  in  Equation  14.  Equation  26  can  be 
written  as: 

S(tk+i)  =  D1z(tk)  +  D2z{tk)  +  5crnd{tk)  [Ds  +  D4]  (27) 
so  that 


D 1  —  i  H-  k*&2,i 

D2  =  ka<&  i52  +  ^^2,2 

D3  =  flk+1  ka$lt 2(4+1  -  r)dr 
fit 


(28) 


£>4  =  £+1  fc$2,2(4+i  -  r)dr 


It  is  now  required  to  evaluate  <f>  and  perform  the  math  re¬ 
quired  in  Equation  28.  By  computing  the  necessary  inverse 
Laplace  transforms,  it  is  found  that 


$1,1 


otf-A  t 


=  sin(uJdAt)  +  cos(uJdAt) 


_V/IZC5 

"  crAt  ,  A  \ 

$1,2  =  ^ sin{udAt ) 
<E>2,i  =  -tvn2^-sin(udAt) 


(29) 


$2,2  =  effAt 


cos(u>dAt) - 7^ — sin(udAt) 

v  !  — C2 


where  odd  and  £  were  defined  after  Equation  16  and  At  = 
tfc+i  —  tk.  Using  Equation  28,  Di,  D2,  D3,  and  £>4  are  found 
to  be 


27 1  =  &ae 


a  At 


C 


U3? 


-;sin(uJdAt)  +  cos(uJdAt) 


0<j  At 


kujn  2 - sin(wdAt) 

y d 

kaeaAt 

D2  =  - sin{uJdAt) 


Ud 


ke' 


a  At 


cos(uOd^t)  — 


c 


-sin(uJdAt) 


K 


V^c2 

ddAt) 

{(jcos((x>dAi)  +  <x>dsm((+>ciAi)}] 

[-^d  -  e<jAt  {<rsm(wdAf)  -  WdCos(wd  Af)}] 


£>3  =  +  e<TA*  {o-sm(wdAt)  -  u)dcos(u)dAt)}] 

Da  =  A  -1-  ^At 


-0  +  e° 


(30) 


The  objective  is  to  find  a  gain  M  that  will  modify 
AScmdcA(tk)  in  such  a  way  that  6(4+1 )  —  $cmdcA  (4). 
Hence,  it  is  desired  to  find  M  such  that 

A  ScrndcA(tk)  +  S(tk)  =  D1z(tk)  +  D2z(tk ) 

+  (D3  +  D4)  [M A5crndCA  (tk)  +  <S(f*)] 

(31) 


Solving  for  M  gives 


M  = 


^ crude A^k)  -^3 

-  27i2?(tfc)  -  D2z(tk) 

(D3  H-  2^4) AScrndc A(tk) 


■  D4)S(tk) 


(32) 


In  this  case,  it  can  be  seen  that  not  only  is  the  actual  actu¬ 
ator  position  needed,  so  are  the  intermediate  variables  z(tk) 
and  z(tk).  In  practice,  an  estimator  (for  example,  a  Kalman 
filter)  would  be  designed  to  provide  z(tk)  and  z(tk).  As 
with  the  first-order  and  simple  second-order  cases,  for  a  sys¬ 
tem  with  more  than  one  actuator,  M  in  Figure  4  would  be 
a  diagonal  matrix  with  entries  along  the  main  diagonal  be¬ 
ing  Mi,  M2,  •  •  •  ,  Mm.  Here,  Mi  would  be  computed  using 
Equation  32  and  ujn . ,  Q  corresponding  to  the  ith  actuator. 

5.  Simulation  Results 

In  this  section,  results  from  a  simulation  of  the  system  dis¬ 
played  in  Figure  4  will  be  shown.  A  rate  and  position  con¬ 
strained  linear  programming  based  control  allocator  will  be 
utilized  in  this  work.  In  this  case,  the  control  allocation  algo¬ 
rithm’s  objective,  referring  to  Figure  3,  is  to  find  Scmd  such 
that 

&des  =  B  Scmd  (33) 

where  B  is  the  control  effectiveness  matrix  and  d des  is  typi¬ 
cally  a  set  of  moment  or  acceleration  commands  for  the  roll, 
pitch,  and  yaw  axes.  Although,  if  feasible,  the  control  allo¬ 
cator  will  be  able  to  find  a  Scmd  such  that  Equation  33  holds, 
the  real  test  is  to  determine  what  happens  after  the  actuator 
dynamics  operate  on  Scmd.  Hence,  the  overall  system  goal  is 
to  achieve  S  such  that 


*des 


=  B  S 


(34) 


—0.4 

0.4 

-0.1 

0.1  ' 

B  = 

-0.1 

-0.1 

-0.6 

-0.6 

-0.1 

0.1 

-0.1 

0.1 

Equation  34  is  the  metric  upon  which  the  quality  of  results 
will  be  judged.  In  this  example,  four  control  effectors  are 
present  and  the  control  effectiveness  matrix  is  fixed  at 


(35) 


where  the  elements  of  B  have  units  of  (rad/ sec2) /deg. 
Since  there  are  more  control  effectors  (4)  than  axes  to  control 
(3),  a  control  mixed  or  allocator  must  be  used.  In  this  work,  a 
flight-tested  linear  programming  based  control  allocation  al¬ 
gorithm  with  rate  and  position  limiting  will  be  used  [2],  [?]. 
It  is  important  to  note,  however,  that  any  control  allocation 
method  could  be  used. 
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In  the  following  simulations,  a  mixture  of  actuator  dynamics 
will  be  used.  Specifically,  the  dynamics  of  each  actuator  are 


<5i  (s)  _  2.5(s+10) 

Semd1(s )  -  s2+7.071s+25 
52(s)  2.50+10) 

Scmd2(s)  -  s2+7.071s+25 
<^3p)  _  49 

(^)  S2+7s+49 

<^4p)  =  5 

ficmd^  0)  s  +  5 


(36) 


which  gives  (  =  uon  =  5,  k  =  2.5  for  the  second-order 
with  zero  cases,  (  =  0.5,  uon  =  7,  k  =  49  for  the  simple 
second-order  case,  and  a  =  10.  Each  actuator  is  rate  and 
position  limited  by  the  following  values 

Smin  =  [  -1.5  -1.5  -1.5  -1.5  ]  (deg) 

Smax  =  [  0.4  1.5  1.5  1.5  ]  (deg) 

SmaxcA  =  [  10  10  10  3  ] 


These  limits  were  selected  so  that  at  least  one  position  and 
one  rate  limit  were  in  effect  at  some  time  during  the  sim¬ 
ulation.  This  was  done  to  show  that  the  method  developed 
in  this  work  is  applicable  when  control  effectors  are  satu¬ 
rated.  As  will  be  shown,  when  compensation  is  used,  actua¬ 
tor  1  becomes  upper  position  limited  and  actuator  4  becomes 
rate  limited  as  the  frequency  of  the  commands  increases. 
The  command  signals,  d^es,  consist  of  chirps  of  magnitude 
0.15,  0.3,  and  0.15  in  the  roll,  pitch,  and  yaw  channels 
and  where  the  frequency  ranged  from  0.5  —  2  Hz  as  a  linear 
function  of  time  over  a  10  sec  time  interval. 


Simulation  runs  were  performed  with  and  without  compen¬ 
sation  for  magnitude  attenuation  due  to  actuator  dynamics. 
Ideal  conditions  are  when  the  actuator  dynamics  can  be  rep¬ 
resented  by  s  ^(g)  =  1.  If  sufficient  control  authority  and 
ideal  conditions  exist,  then  the  control  system  would  achieve 
d des  =  B5,  which  is  the  best  possible  performance.  Fig¬ 
ures  5  and  6  show  d des  and  B<*>  without  and  with  the  mag¬ 
nitude  compensation  described  in  Equations  7,  19,  and  32 
applied.  Clearly,  when  the  magnitude  compensation  is  not 
used  (Figure  5),  d des  ^  B<5,  and  a  large  error  exists  between 
these  two  quantities.  When  magnitude  compensation  is  used 
(Figure  6),  however,  d^es  =  Bd  and  near  ideal  performance 
is  achieved. 


Figures  7  and  8  show  the  control  effector  commands, 
<5 cmdc a’  and  actual  deflections,  d ,  when  magnitude  compen¬ 
sation  is  not  utilized.  As  expected,  the  actual  deflections  are 
considerably  different  than  the  commanded  deflections  and 
this  directly  translates  to  the  deviations  between  d des  and  Bd 
as  shown  in  Figure  5.  On  the  other  hand,  Figures  9  and  10 
display  the  control  effector  commands,  ScrndCA ,  and  actual 
deflections,  S ,  when  magnitude  compensation  is  used.  In  this 
case,  the  actual  deflections  are  nearly  equal  to  the  commands 
and  as  a  result,  d des  —  B<5.  Notice  that  control  effector  1  is 
upper  position  limited  and  control  effector  4  is  rate  limited 
for  the  last  few  seconds  of  the  simulation  run.  Thus,  using 


Figure  5.  Commanded  (d^es)  and  simulated  (B£)  angular 
accelerations  -  Compensation  OFF. 


Time  (sec) 


Figure  6.  Commanded  (d^es)  and  simulated  (B£)  angular 
accelerations  (jj$)  -  Compensation  ON. 


the  simple  gain  adjustment  described  in  Equation  19  results  in 
the  actual  control  deflections  being  equal  to  the  commanded 
control  deflections.  It  is  apparent  that  adjusting  the  control 
effector  command  increments  can  help  to  mitigate  adverse  in¬ 
teractions  between  discrete  time  implementations  of  control 
allocation  algorithms  and  actuator  dynamics. 

6.  Conclusions 

Interactions  between  constrained  control  allocation  algo¬ 
rithms  and  the  dynamics  of  actuators  can  result  in  degraded 
performance  if  not  carefully  implemented.  It  was  shown  that 
the  control  effector  commands,  from  a  control  allocation  al¬ 
gorithm,  are  attenuated  by  actuator  dynamics.  One  method, 
which  can  be  used  to  extract  maximum  performance  from 
such  a  system,  is  to  modify  the  control  effector  command 
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Figure  7.  Control  Effector  1  and  2  Positions  - 
Compensation  OFF. 


Figure  9.  Control  Effector  1  and  2  Positions  - 
Compensation  ON. 


Figure  8.  Control  Effector  3  and  4  Positions  - 
Compensation  OFF. 


increments  from  the  control  allocator.  The  gains  used  to 
modify  the  commands  were  computed  for  first-order,  simple 
second-order,  and  second-order  with  a  zero  actuator  dynam¬ 
ics.  Simulation  results  show  that  significant  improvements 
can  be  achieved  by  using  this  method  to  post-process  the 
commands  from  the  control  allocation  algorithm.  Benefits 
of  this  method  are  that  the  control  allocation  algorithm  need 
not  be  modified,  it  requires  minimal  additional  computations, 
and  it  is  valid  for  saturated  and  unsaturated  control  effectors. 
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